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^ ■ Abstract 



A microscopic theory of the electromagnetic radiation emitted by a highly excited 
f-^. ; nucleus is developed on the basis of the Landau theory of a Fermi liquid. Closed 

■ formulae are obtained for the mean radiative width and its mean square fluctuation 

from level to level. The temperatures of many nuclei are found from the observed 
widths. The relaxation time is estimated from the experimental data on the radiative- 
width fluctuations. The regions of applicability of the various types of relations between 
the relaxation time and the lifetime of the compound nucleus, as well as the relevant 
physical consequences, are discussed. 

"o . 

§ ; 1 Introduction 

That there is a gradual unification of the mechanisms underlying the emission of electromag- 
netic radiation by a not too light nucleus as the excitation energy of the nucleus increases 
is practically beyond question. Roughly speaking, if an infinite number of the levels of the 
nucleus as a whole lie below the initial excitation energy, then the system will itself find 
and prefer the easiest effective way of emitting 7 quanta. In particular, there is below the 
initial excited level an abundance of levels 7 transitions to which satisfy the most favorable 
"selection rules" , so that from this point of view the process is, in the limit under consider- 
ation, practically one with an infinite number of channels. Analysis shows (see [1]) that the 
electric-dipole radiation due to collisions between the (proton and neutron) quasiparticles 
and the "wall" of the nucleus predominates. The ideas of the Landau theory of the Fermi 
fluid [2,3] allow us to compute in closed form the radiative width T 7 and its fluctuations 
from level to level. 

In a sense, spherical nuclei are rather exotic objects: the application to them of the 
Fermi-fluid concepts requires in each specific case certain precautions. The fact that the 
spherical configuration is stable is in itself an indication of the essential role played by the 
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"residual interaction" between the quasiparticles, an interaction which blurs the Fermi level: 
it can be shown that in the scheme without interaction the sphere is absolutely unstable 
(see [4]). Furthermore, analysis of the data on the shell and magic oscillations in the masses 
of spherical nuclei allows the establishment of the macroscopically ordered structure that 
this residual interaction possesses in the space of the values of the orbital momentum / 
of the individual quasiparticles [5]. (We shall again touch upon spherical nuclei when we 
compare below the theoretical results with the experimental data.) Nonspherical nuclei, on 
the other hand, are easy to investigate, it being apparently necessary to regard their shape as 
a perfectly natural consequence of the properties of the "normal" , disordered nuclear phase, 
in which the quasiparticles situated near the limit of the Fermi distribution move, in the 
main, independently of each other. 

However, for the theory of radiative widths expounded below the "shape effects" , as such, 
are of no independent importance, and a special allowance for them is not necessary. Indeed, 
the equilibrium deformation a of a nonspherical nucleus is equal in order of magnitude to 
pj 1 (i.e., a ~ P/ 1 )) where 

p f = k f R^>l (1) 

(kf is the limiting momentum of the quasiparticle distribution and R is the radius of the nu- 
cleus) is an important dimensionless parameter that arises in the most diverse investigations 
that have as their aim the treatment of the nucleus as a macroscopic body. In view of the 
scalar nature of the quantity T 7 to be computed, only the squared deformation can enter, so 
that the relative magnitude a 2 ~ p"j <C 1 of the corresponding corrections is negligibly small 
and falls outside the limits of accuracy of the theory. In other words, only the possibility 
of considering the quasiparticles individually is important for what follows, it still being 
possible in actually occurring deformations to treat the geometry of the motion of each of 
the quasiparticles as spherically symmetric. 



2 The radiative width of a highly excited nuclei 

We shall derive the expressions, referred to one quasiparticle, on the basis of the corre- 
spondence principle. As applied to radiation processes, this principle asserts total analogy 
between the formulas of the classical and quantum theories (see, for example [6]). The clas- 
sical intensity / (i.e., the energy emitted per unit time) needs only to be divided by e = Jtlu 
to be converted into the quantity of real interest — the probability of emission of individual 7 
quantum. The only remaining difference consists in the following: the spectral component of 
the multipole moment, which varies according to a classical law, should, generally speaking, 
be replaced by the corresponding matrix element of its operator. They, however, coincide in 
the quasiclassical limit (see (1)) [7]. Consequently, we can speak of a quasiparticle trajectory: 
it is in this case a straight line joining two points on the surface of the nucleus. 

The basic formula of the classical theory of the electric dipole radiation has the form 

(2) 

where is the charge of the radiating particle (the radiating quasiparticle, in the general 
case; see below) and r is its acceleration vector. Before proceeding to the spectral decompo- 
sition of /, let us note that in the thermal-equilibrium state the quasiparticle motion inside 
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the nucleus does not change its qualitative character in time. Therefore, let us formally carry 
out the Fourier expansion over an arbitrary, but sufficiently long interval of time t and then 
take the limit as t — > oo (see [8]): 
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In view of the independence of the different chords traced by the quasiparticles in their 
wall-to-wall motion, we have 
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where t\ = l/vf is the time it takes to travel from one end of a chord to the other, I = 
2 a/-R 2 — p 2 is the length of the chord, and n is the number of chords. Here we have taken 
into account the fact that in the region of the radiation energy spectrum of interest to us 



out! < 1 



(5) 



(we shall return to the criterion (5) later). Furthermore, here and below the quasiparticle 
velocity v is replaced everywhere by its limiting value Vf. The point is that because of the 
Pauli principle only those "elementary emitters" (i.e., quasiparticles) that are situated in 
the immediate neighborhood of the Fermi level play a role (see below). 

The distribution of the chords over the impact parameters p is easily found from consid- 
erations of isotropy and homogeneity of nuclear matter: 



3 rR 

w{p)dp= ^V R2 ~ P 2 pdp; J w(p)dp=l. 



(6) 



Averaging, in accordance with (4) and (6), the square of the integral over the radius vector, 
we also express the number n of quasiparticle — nuclear wall collision events in terms of the 
physical time t: 
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Taking into consideration the relations (2)-(7) and the relevant considerations, we obtain 



f(e)de = 
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(8) 



This expression gives the probability per unit time of emission by one quasiparticle of a 7 
quantum in the interval de of its energy values. 

According to the theory of the Fermi fluid [2-3], the mean occupation numbers of the 
individual quantum states are given by the standard Fermi distribution 
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where e' is the quasiparticle energy measured relative to the chemical potential and T is 
the temperature. On the other hand, the number of actual single-quasiparticle states in the 
volume V = (4/3)7i\R 3 is equal to 

— \f tP 1 4 R 3 o 2 

(10) 

(pf = hkf is the limiting momentum in standard units), where allowance has been made for 
the additional spin doubling. In fact, even in the quasiclassical limiting case, (1) remains an 
important quantum effect due to the identity, the indistinguishability of identical fermions 
[7] : the above-described classical picture of the process is actually realizable only in the case 
of radiative transitions that are compatible with the Pauli principle. Therefore, the product 
of the expressions (8), (9), and (10) should be supplemented by the factor 

1 — n(e' — e), 

which determines the fraction of the transitions admissible by this principle. Then integration 
over the energies e' of the radiating quasiparticles reduces to 

nW[l-n(e'-e)]de'=^- r[ . (11) 

Integration over the boson energies yields 

poo _4 
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where ( is the Riemann zeta function. Finally, the 7-quantum emission probability per unit 
time, appropriately summed over the entire set of quasiparticles of the same sort will be 
given by 

W - —C(5)^R 6 T 5 (13) 

where m* — Pf/vf is the effective mass of the quasiparticle. 

Above, as the coordinate origin convenient for the calculations, we used the geometri- 
cal center of the nucleus. However, the role of the total charge Ze of the whole system in 
processes induced by the oscillations of the radius vector of the individual nucleons (quasi- 
particles) is well known. Because of recoil, even the electrically neutral quasiparticles (i.e., 
the neutron quasiparticles) will appear to emit radiation during their motion relative to the 
center of the nucleus. The corresponding, well-known, "charge-renormalization" formulae 
have the form 

4 = (l - |) e, e N d = -|e (14) 

(the "effective charges," (14), of the two components are correct only for processes induced 
by the oscillations of the electric dipole moment of the nuclear system (see, for example, 
[7])). Summing, with allowance for (14), the expressions (13) or the proton and neutron 
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components of the nuclear matter, and multiplying them by H in order to convert them into 
the energy widths of interest, we finally obtain 



1024 eV[ Z( ZY\ 
^ = 105^ C(5) W V~ 2 A \ l ~A)\ RT (15) 

(notice that the numerical factor (1024/1057r 2 )£(5) = 1 is very close to unity). The law 
T 7 oc T 5 was given in the preceding paper (see [1], formula (4)), where it was motivated by 
semi-phenomenological considerations. 

To what extent can the result (15) be identified with the radiative widths of the individual 
resonance levels of a compound nucleus that is excited, say, in a reaction involving slow- 
neutron capture? It follows from its derivation that the formula (15) corresponds to a state 
in which at the temperature T the quasiparticles of the nuclear Fermi liquid are in thermal 
equilibrium with each other. On the other hand, the width T 7 of a specific level can, reasoning 
abstractly, be conceived to have been computed from some very complicated, unknown (to 
us) wave function of the corresponding state of the nucleus as a whole. According to the 
fundamental principles of statistics, the two approaches lead to results that coincide to within 
the values of the fluctuations (see, for example, [3]). 

Let us now rewrite the condition (5) of applicability of the theory in a more concrete 
form. Owing to the thermal nature of the radiation, the inequality (5) is equivalent to the 
following inequality: 

T<e , (16) 

where 

e = hv f /2R ~ 5MeV (17) 

is the characteristic energy corresponding to the reciprocal of the time it takes a quasiparticle 
to cross the nucleus along a diameter. 

It is worth noting that in the opposite limiting case 

T>e , (18) 

because of the oscillations of the exponent e %u)t along the chord traced by the quasiparticle 
(see (3) and (4)), the energy distribution of the 7-quantum emission probability acquires the 
form of the well-known Planck black-body radiation spectrum [3]. The radiation width 
would, accordingly, become proportional to the cube of the temperature in the case of 
a sufficiently strict fulfilment of the condition (18). However, this "black-body radiation 
limit" defined by (18) is, in practice, hardly attainable in nuclear physics. At least the 
temperature of the compound nucleus should not exceed the nucleon binding energy, which 
is ~ 8 MeV — otherwise the neutrons would fly out of the nucleus "instantly," escaping the 
thermal-equilibrium establishment phase 1 . 

1 Notice that the radiation due to collisions between the quasiparticles would then become dominant only 
at T > Jpje^. 
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3 The mean-square fluctuation in the radiation width. 
The role of the relaxation time 



The direct, quantum- mechanical computation of the characteristics of the individual states of 
the nucleus is inexpedient and practically impossible. Furthermore, as applied to macroscopic 
bodies (see the criterion (1)), this, as a rule, borders on the theoretical impossibility [3]. 
Therefore, we are obliged here to treat the state of the occupation of the individual quantum 
states of the quasiparticles of the Fermi liquid as a randomly varying function of the time. We 
shall calculate the instantaneous "emissive power" T 7 of the nucleus in a manner completely 
similar to the computations of the preceding section. The "one-component" variant of the 
corresponding formula can be represented in the form 
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(19) 



Here we have, for simplicity and convenience , written the discrete sum over the fermion 
states. In case of need the transition to integration can easily be accomplished with the aid 
of (10). 

The "instantaneous" , physically realizable values 



0, 1 



(20) 



of the fermion occupation numbers differ from the mean occupation numbers (9). This 
circumstance is the obvious cause of fluctuations in Fermi systems. It is convenient to 
consider them with the aid of the simple relation (see [3]) 
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(21) 



Let us find the mean-square fluctuation of the expression (19) — the number of summations 
and integrations doubles upon squaring. One summation over the quasiparticle states is 
trivial owing to the presence of the 5 symbol on the right-hand side of (21); the subsequent 
integration is elementary, although somewhat tedious. Adding, in accordance with (14), 
the squares of the fluctuations in the proton and neutron components and introducing the 
dimensionless variables = £ i,2/T in place of the 7-quantum energies, we obtain 
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The details of the integration over the boson energies are given in the Appendix. The final 
result has the form 



j 848 8 
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(23) 



n=l ' ' k=n+l 

Notice that the term with the double sum is about half percent of the value of the integral, 
so that in practice we can restrict ourselves to the consideration of only the first term on the 
right-hand side of (23). 
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The problems pertaining to the fluctuations are relatively subtle and require a more 
careful physical treatment. In particular, there is no reason to equate (Ar 7 ) 2 to the mean 
square (Ar 7 ) 2 of the actually observable, physical fluctuation in the radiative widths of 
many close resonance levels. This becomes especially apparent when we consider the most 
important and interesting case in which thermal equilibrium in the nucleus is established 
long before the "decay" of the nucleus: 

IY//i«l. (24) 

Here r is the relaxation time (see below) and T is the total width of the initial state of the 
nucleus. Taking into account the fact that this quasistationary state decays according to the 
law exp(— rt/H), we express the number v of emitted quanta and its fluctuation in terms of 
the instantaneous emissive power T 7 : 

1 f°° ~ 

v = ^-\ r 7 (t) e - r '/^t, 

* Jo 

(Az/) 2 = ^ J J™ Af 7 (t)Af 7 (t') exp [-r(t + t')/h] dtdt', ( 25 ) 

Af 7 (t) = f 7 (t) - r 7 . 

Further, it is convenient to introduce the notation t' — t + r. According to the thermo- 
dynamic theory of non-equilibrium processes (and of the corresponding fluctuations in the 
thermal-equilibrium state; see, for example, [3]), the mean value of the time correlation of 
the fluctuations is given by the relation 



Ar 7 (t)Ar 7 (t + r) = (Ar 7 ) 2 exp(-|r|/f), (26) 

where f is the relaxation time. With allowance for (24), the substitution of (26) into (25) 
yields 



Let us now consider the ensemble of the large number of close levels of a compound 
nucleus with radiative width V: owing to the fact that the levels decay according to the 
single law exp(— Tt/h), the equilibrium in the ensemble (the equipopulation of the levels) is 
not destroyed in time. The number of 7 quanta 

v = I r T.e-^dt = ^ 

frJo 7 r 

has been pre- averaged over a group consisting of many levels with practically the same T 7 . 
In the final averaging of the square of the fluctuations (Az/) 2 over the entire ensemble of the 
groups differing in their radiative widths T 7 , each group is taken into account with a weight 
proportional to the number of levels in it: 



(Az/) 2 = (Ar 7 ) 2 /r". (28) 
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Equating the right-hand sides of the formulas (27) and (28), we finally obtain 



(Ar,)^ = (Tr/h)(AT,y. (29) 

(we shall no longer write the averaging sign over the relaxation time r). A striking feature of 
the relation (29) consists in the following: It turns out that the fluctuations in the probability 
of decay of the compound nucleus via the radiative channel depend on the total decay 
probability T, including all the generally possible decay channels. The physical meaning 
of the formula (29) is simple: In the time picture the deviation of AT y (t), the emissive 
power, from its mean value has time to average out to some extent provided the decaying 
exponential function varies sufficiently slowly (see the criterion (24)). The small factor Tr/h 
on the right-hand side of (29) is precisely the quantity that determines the fraction of the 
physical, actually observable effect that remains after such a partial averaging. 

4 Comparison with experiment 

With the aid of the formula (15) we determined the temperatures of compound nuclei from 
the observed radiative widths of their resonance levels [9,10]. The results of such an analysis 
for two well-known regions of nonspherical nuclei are given in the tables. We assumed in the 
computations that 

R = 1.2 • 10~ 13 A 1/3 cm (30) 

and m* = m n , where m n is the mass of the free nucleon. It is noteworthy that the tempera- 
ture in the case of the actinide nuclei turns out consistently to be ~ 100 keV lower than the 
characteristic temperature for the lanthanide region. This may be due to both the decrease 
of the neutron attachment energy toward the end of the Mendeleev periodic table and the 
difference in the atomic weight A. A similar temperature decrease apparently occurs within 
the nonspherical-lanthanide region. 

Spherical nuclei possess a number of unique features that must be taken into consideration 
(see the Introduction). However, the question of the applicability to them of the formula 
(15) is at present difficult to answer categorically. Indeed, nuclei of this sort apparently 
undergo a phase transition to the "normal" , nonspherical state at temperatures 

T ~ Ae', 

where Ae' is some characteristic width of the diffuse zone of the Fermi distribution, a zone 
which owes its existence to the residual interaction. Meanwhile, the spectrum of the emitted 
quanta (it is given by the integrand on the left-hand side of (12)) is such that the energy 
averaged over it is equal to 

e = 5T 

(see also [1], formulae (5) and (6)). Thus, many of the radiative transitions can, roughly 
speaking, elude that region of the statistical distribution of the quasiparticles where the 
distribution differs significantly from (9). Therefore, the attempts to apply the formula (15) 
also to spherical nuclei, though not rigorous, is nevertheless of some interest. It is natural 
to suppose that spherical nuclei have higher temperatures (and, consequently, relatively 
low entropies; see also [1]). Comparison with the data on the radiative widths apparently 



8 



Compound nucleus 


E max , MeV 


T 7 , MeV 


T, MeV 


Nonspherical lanthanides 




62 Sm^ 8 


8.14 


52 


0.42 


Q™150 

62om 88 


7.98 


64 


0.44 


C m 153 


5.89 


71 


0.45 


"Rn 152 

63^U 89 


6.29 


89 


0.47 


6 3 Eu 9 f 4 


6.39 


102 


0.48 


p J 156 


8.53 


110 


0.48 


64Gdgg 7 


6.35 


110 


0.48 


PJ158 


7.93 


89 


0.46 


pjl59 


6.03 


105 


0.48 


TV. 160 
65 1 D95 


6.40 


90 


0.47 


66L>y 9 6 


8.20 


122 


0.49 


66Dyg7 3 


6.25 


175 


0.52 


66^798 


7.66 


103 


0.47 


66Dygg 2 


5.64 


166 


0.51 


66DY96 2 


8.20 


122 


0.49 


67nOgg 


6.33 


91 


0.45 


68^ r 99 


6.44 


97 


0.46 


rp^l68 
68^100 


7.77 


96 


0.46 


Trn 170 

69-Lm 101 


6.38 


86 


0.44 


70 I D 102 


8.14 


74 


0.43 


70 i D 104 


7.44 


79 


0.43 


72Hfioe 


7.62 


64 


0.41 


T 182 


6.06 


54 


0.39 


wl83 
74 W 109 


6.19 


58 


0.40 


74VV 110 


7.42 


74 


0.42 


AI7185 

74 w 111 


5.75 


64 


0.41 


74 W 113 


5.46 


62 


0.40 


D„188 

75 ite 1131 


6.24 


55 


0.39 




5.73 


55 


0.39 
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Compound nucleus 


E max , MeV 


r 7 , MeV 


T, MeV 


Nonspherical actinides 




rpVi 233 
90-Ln 143 


4.96 


21 


0.30 


Pq232 


5.52 


44 


0.33 


T~) 234 
91^ a 143 


5.12 


48 


0.35 


TT234 
92 U 142 


6.78 


40 


0.34 


T t235 
92 U 143 


5.27 


25 


0.31 


T T236 
92 U 144 


6.47 


40 


0.33 


TT237 
92 U 145 


5.30 


29 


0.31 


92 U 147 


4.78 


23 


0.30 


93 Np?I 


5.43 


34 


0.32 


94ru 146 


6.46 


40 


0.33 


ID 241 

94 ru 147 


5.41 


31 


0.32 


p.. 242 
94^U 148 


6.22 


37 


0.33 


Am 242 

95 Am 147 


5.48 


42 


0.33 


Am 244 
95 Am l49 


5.29 


50 


0.35 


96^ m 148 


6.72 


37 


0.33 


96Cm^ 4 g 


5.70 


39 


0.33 


96Cm^ 51 


5.21 


35 


0.32 



corroborates this trend. For example, for the compound nucleus 79 Au^g we obtain T = 0.45 
MeV, in the case of 80 Hg 2 ° 2 we have T = 0.56 MeV, and, finally, T = 0.62 MeV for 8 iTl 2 ^. 

The experimental study of radiative-width fluctuations became possible only recently 
as a result of an increase in the accuracy of their measurement, and comparison of the 
theoretical formulas with experiment meets for the present with certain practical difficulties. 
Let us discuss three specific nuclei, for which a selection of resonance levels with accurately 
measured radiative widths nevertheless allowed the estimation of the relaxation time r from 
the formula (29) (see also (22) and (23)). Data on two gadolinium isotopes are given in [11]; 
we took into consideration only the levels for which the error in the radiative width is less 
than 10" 2 eV. In the case of 64 Gd^ 6 (10 levels) f 7 = 0.11 eV, [(Ar\) 2 ] 1/2 = 0.016 eV, and 
h/r = 26 eV. For 64 Gd£f (11 levels) we have f 7 = 0.089 eV, [(AT^) 2 ] 1 / 2 = 0.0087 eV, and 
H/t — 52 eV. Let us also give the results of a similar analysis of the data on holmium [12]: 
87 Ho^ 6 (21 levels) f 7 = 0.091 eV, [(Ar 7 ) 2 ] 1 / 2 = 0.0099 eV, and H/t = 39 eV. Thus, as far 
as we can judge, h/r ~ 50 eV and r ~ 10~ 17 sec, which is a remarkably long time on the 
nuclear scale. We must, however, not forget that the longest of the relaxation times r of the 
system enters into the thermodynamic theory (see formula (26)). The "partial equilibrium" 
at each moment of time was understood to have been established over the significantly 
shorter relaxation times 2 . 

2 From the data of the recent paper [13] we find that H/t <~ 20 eV for the compound nuclei Th 233 and 
U 239 . 
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5 Discussion. Is there enough time for the establish- 
ment of thermal equilibrium in a nucleus? 

The question of relaxation in nuclear matter is of considerable interest. Thus far, as far as 
we know, it has not been possible to estimate the characteristic time of this process on the 
basis of any direct analysis of the experimental data. Therefore, the observed fluctuations in 
the radiative widths can be a valuable source of such information, and even the preliminary, 
tentative figures (r ~ 10 -17 sec; see the preceding section) need to be discussed. For the 
above-mentioned particular cases of resonance excitation by neutrons of energy < 1 keV, the 
condition (24) was satisfied with three orders of magnitude to spare — in other words, total 
thermal equilibrium was attained in the nucleus. The situation can, however, change when 
we go over to higher kinetic energies of the bombarding particles (see below). 

It would, apparently, be somewhat naive to regard the time t\ = H/eq ~ 10~ 22 sec of 
transit of a quasiparticle through the nucleus as an estimate for the relaxation time. Indeed, 
for example, the system is not in the least drawn nearer to the state of thermal equilibrium by 
a coherent, reversible, purely elastic act of collision between a quasiparticle and the surface 
of the nucleus 3 . On the other hand, inter-quasiparticle collisions appear to be quite an 
effective relaxation mechanism; it may well turn out to be the dominant mechanism. 

For obvious reasons, the thermal-radiation data used in the present paper actually per- 
tained to an extremely narrow part of the energy spectrum of the system. Let us now 
qualitatively consider how the real conditions under which the relaxation process proceeds 
in the nucleus should change upon further increase in the excitation energy of the nucleus. 
The neutron width T n increases first in proportion to the square root of the distance from 
the neutron-detachment threshold; then there we come into the stage of much more rapid 
exponential growth. It is well known that owing to this phenomenon the resonance levels 
merge, forming a continuous spectrum. But then whether the lifetime h/Y of the compound 
nucleus will be long enough for the establishment of total thermal equilibrium in the nucleus 
may become doubtful, beginning from 4 

r = T n > 100 eV. 

In this connection, it is desirable to try critically reinterpret the method, based on the 
neutron "evaporation" process, for determining nuclear temperatures. It is difficult for the 
present to judge how the fact that the state of the neutron-emitting nucleus is not a totally 
equilibrium state will influence such an analysis. Not much doubt has thus far been expressed 
about the evaporation temperatures probably because their order of magnitude is quite 
plausible (and, in so far as we can judge, indeed correct). However, as the experimental 
investigation of the reactions (n, n') goes on, attention will have to be paid not only to the 
absolute figures, but also to the behavior of the relevant quantities. Of special interest, 

3 Besides, a characteristic time ~ t\ is quite capable of playing an important role at the earliest stage 
of the development of the nuclear reaction. Here, however, we are discussing only the late, final phase of 
the thermal-equilibrium establishment process — see the end of the preceding section. In particular, it is 
extremely doubtful that there will, at such times, remain reasonable physical criteria for distinguishing the 
initial bombarding particle (or the corresponding quasiparticle). 

4 It is possible that the so-called Erickson fluctuations in the nucleon-nucleus interaction cross sections 
[14] are also due to this circumstance. 
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in particular, is the case when the temperature of the compound nucleus as a function of 
its excitation energy is an almost horizontal, non-monotonic in detail, and often simply 
a decreasing function. Unfortunately, the authors of the corresponding publications give 
this remarkable circumstance comparatively little consideration (see, for example [15]). We 
could have attempted to interpret the decrease of the temperature with increasing excitation 
energy as some giant random fluctuation, but it would have been difficult to conceive it as 
a phenomenon that would occur with any degree of consistency. On the average, however, 
negative specific heat is impossible for the nucleus. A state with negative specific heat is 
totally unstable, and cannot be realized in nature [3]. 

We express our thanks to A.I. Baz', V.K. Voitovetskii, LI. Gurevich, A.G. Zelenkov, L.P. 
Kudrin, V.M. Kulakov, A.A. Ogloblin, I.M. Pavlichenkov, N.M. Polievktov-Nikoladze, and 
K.A. Ter-Martirosyan for a discussion of the results of the paper. 



Appendix 

Let us expound in some detail the integration over the boson energies in the formula (22). 
The terms of the integral J are identical in form, but individually each of them contains a 
pole at X\ — x 2 - Therefore, it is sufficient to evaluate any of these integrals in the principal 
value sense: 

J = 2 [°° dx 1 -x 4 1 e Xl coth^P [°° ^ -. (A.l) 

Jo 2 Jo (e* 1 -e**)(e*i - e~ x *) 

Let us transform the inner integral with the aid of the substitution y = e~ X2 , representing it 
as a derivative with respect to some parameter: 



p r x\dx 2 = £ _ x p r 1 (Inyfdy 

J Q ( e *i -e* 2 )( e *i - e ~**) J Q 



o (y-e x )(y-e x ) 
ov Jo {y-e x )(y-e x ) 



(A.2) 



(x — Xi). In the decomposition 

^-^-^ ^W' 1 -^'"' (A - 3) 

in which the integrand is expressed in partial fractions, the integrals 



o y~e x Jo y-e- 



can conveniently be expressed as series. For this purpose, let us represent the fractions by 
the corresponding geometric progressions, and let us also take into consideration the formula 



OO j 

= 7T ( 



n=—oo 



" COt 7TU. 

n + v 
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We then have 

— TLX 1 — 7132 

n=l n=l 

Taking the limit in accordance with (A. 2), we obtain 

J = 192 i„ \3-T« "go 5 }^ + 6 g^X 5S?^ d «- < A ' 6 > 

It is easy to see that the integral standing under the summation sign 

^ f'MV^/ p f' _ _ (A . 7) 



o sinh 2 (e/2) J (y-1) 2 J ( y - e*)(y - e-*) 

is essentially of the same type as (A. 2), except that it is differentiated once more with 
respect to v and that the additional passage to the limit x — > will also be necessary. The 
subsequent simple computations, besides the passages to the limit, also include convenient 
re-designations of the indices of the double summation. As a result, we obtain the formula 
(27). 
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